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ABSTRACT 

High-resolution HST imaging of the compact planetary nebula NGC 6644 has 
revealed two pairs of bipolar lobes and a central ring lying close to the plane of 
the sky. From mid-infrared imaging obtained with the Gemini Telescope, we 
have found a dust torus which is oriented nearly perpendicular to one pair of the 
lobes. We suggest that NGC 6644 is a multipolar nebula and have constructed 
a 3-D model which allows the visualization of the object from different lines of 
sight. These results suggest that NGC 6644 may have similar intrinsic structures 
as other multipolar nebulae and the phenomenon of multipolar nebulosity may 
be more common than previously believed. 

Subject headings: ISM: general — planetary nebulae: individual (NGC 6644) 
— stars: AGB 



-2- 



Introduction 



Although planetary nebulae (PNs) are usually assumed to have a simple spherical ge- 
ometry, the actual morphologies of PNs have diverse forms. When deeper, narrow-band 
observations are taken, PNs often reveal outer structures such as lobes and haloes. While 
it has been recognized for some time that there are many examples of planetary nebulae 
with prominent bipolar structures (e.g., NGC 2346, NGC 6302), an increasing number of 
nebulae with multiple lobe structures have been found (e.g., NGC 2440, NGC 6072). The 
origin of such structures and the physical mechanisms leading to their formation are yet to 
be understood. 

Since PNs are made up of ionized, molecular, and dust components, a comprehensive 
understanding of the dynamical processes that create the observed morphological structures 
can only be obtained with multiwavelength observations. With modern 8-m class telescopes 
equipped with high-resolution mid-infrared imaging capabilities, the distribution of the dust 
component can be compared to the distribution of the ionized gas component with com- 
parable resolution. In particular, the infrared observations will allow us to search for the 
existence of a high-density torus, and determine what role the torus play in the collimation 
of observed optical outflows. 

The P N NGC 6644 (PNG 008.3-07.3, Hen 2-408, VV 188, IRAS 18295-2510) was dis- 



covered by iHubbld ( I192ll ) . Based on early ground-based imaging, this compact nebula has 



been classified as an elliptical PN by iStanghellini et al.l ( 119931 ). It has a la rge heliocentric 



1999 


). However, 


Gonv et al. 


(2004) 



on their spectroscopic study. Although NGC 6644 has a relatively high surface brightness, 
it has received very little attention in the literature. There are only three entries in ADS on 
the object after a search of keyw ords in the abs t ract, and the papers are mostly related to 



abundance studies of the nebula. lAller fc Keyesl (119881 ) proposed that the object was made 



from a much less metal-rich mixture than the Sun. 

In this paper, we present results of our morphological study for NGC 6644 based on 
optical and infrared images taken with HST and Gemini Telescope, respectively. These 
images reveal for the first time unique morphologies and properties of the system. The 
observations and data reductions are described in § 2. In § 3, we present the results of 
imaging and spectroscopy in the visible and mid-infrared for this nebula. A discussion of the 
results, and a comparison of the observations with a three-dimensional model are presented 
in § 4. Finally, a conclusion is given in § 5. 
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2. Observations and Data Reduction 

2.1. HST Narrow-band Imaging 

Our HST data on NGC 6644 were retrieved from the Space Telescope Science Archive. 
The images of NGC 6644 were obtained under program 8345 (PI: R. Sahai) using the Wide 
Field Planetary Camera 2 (WFPC 2) on HST. The object was observed with the Planetary 
Camera (PC) on February 26, 2000, which provides a 36."8 x 36."8 field of view (FOV) 
at a spatial resolution of 0.045" pixel -1 . The actual observations were made with different 
exposure times (from 20 s to 400 s) to allow for the imaging of both the bright central region 
and the faint outer parts. The data were processed through the HST pipeline calibration. 
Standard bias subtraction and flat-field correction were performed. Data were taken in two- 
step dithered positions to enhance spatial sampling and cosmic rays removal by using the 
task crrej in the STSDAS package of IRAF. The processed F656N (Ha) observation (X p = 
6564 A, AA = 22 A) with a total exposure time of 1120 s is shown in Fig. [TJ 



2.2. Gemini T-ReCS observations 

The imaging observations were made with the Thermal-Region Camera Spectrograph 
(T-ReCS) on the Gemini South Telescope under program GS-2007A-DD-9 on 2007 October 
10. The field of view of the detector is 28."8 x 21."6 with a pixel scale of 0."09. The 
background emission from the sky and telescope were removed with chopping and nodding 
during the observations. NGC 6644 was imaged with three narrow/medium band filters: 
[S iv] (A c = 10.52 /mi, A A = 0.17 //m), PAH-2 (A c = 11.30 /an, A A = 0.61 /zm), Si-5 
(A c = 11.66 /mi, A A = 1.13 /mi), and the broad Q a band filter (A c = 18.30 /mi, A A = 
1.38 /mi). The total on-source exposure times for the [S iv], PAH-2, Si-5, and Q a filters 
are 455.1, 461.5, 238.9, and 238.9 s, respectively. The measurements were made under good 
(20 percentile) sky conditions, and the mid-IR images are diffraction limited. The measured 
FWHM values for the standard stars are 0.375, 0.382, 0.427, and 0.552 "for the [S iv], PAH-2, 
Si-5, and Q a filters respectively. Photometric standard stars were observed for t he purpose of 



flux ca libration using the spectrophotometric fluxes for the standard stars from ICohen et al. 



( 119991 ). The journal of observations is summarized in Table 1. 
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2.3. Spitzer IRS Spectrum 



The mid-infrared spectra of NGC 6644 were obtained by the IRS (IHouck et al.l 120041 ) 
through the observation program 3633 (PI: M. Bobrowsky) with Astronomical Observation 
Request (AOR) key of 11334400. The measurements were observed using the Short-Low 
(SL) module (5.2 /iin - 14.5 /mi), Short High (SH) module (9.9 /im - 19.6 /mi), and the 
Long High (LH) module (18.7 /im - 37.2 /im) with spectral dispersions of R ~ 600. The 
diaphragm sizes are 3"6 x 57", 4"7 x 1T/3, and ll'.'l x 22'.'3 in SL, SH and LH modules 
respectively. The total integration time of IRS observation is 347 s. The IRS measurement 
of NGC 6644 was not positioned at the central star but on the outer lobe. An overlay of the 
IRS apertures on the optical image of the nebula is shown in Fig. [2J 

Data were reduced starting with basic calibrated data (BCD) from the Spitzer Science 
Center's pipeline version sl8.7 and were run throu gh the IRSCLEAN program to remove 



rogue pixels. Next the SMART analysis package (IHigdon et al.l I2004J ) was used to extract 
the spectra. A final spectrum was produced using the combined IRS observations to improve 
the signal-to-noise ratio (S/N). Since the IRS spectrum in the short wavelength range is taken 
from the outside of the nebula, some scaling is needed for the shorter wavelength observations. 
In Fig. [3j we scale the IRS SH and SL observations by factors of 1.29 and 1.625 respectively 
and are able to obtain a smooth spectrum. When we compare these adjusted flux levels 
to the photometric measurements obtained from the 4 Gemini T-ReCS filters, we find good 
agreement (Fig. [3]). Further confidence of this adjustment is found by the agreement between 
the spectra and the IRAS 12 and 25 /im fluxes. Also plotted in Fig. [3] are the 4 Gemini filter 
transmission functions. We can see that the [S iv] and PAH-2 filters have most of their fluxes 
contributed by the [S iv] line and the 11.3 /im aromatic infrared band (AIB) respectively. 
The Si-5 filter mainly takes in fluxes from the 12 /mi plateau emission feature and the Q a 
filter measures the general dust continuum. Although the Q a filter covers the [S ill] line, we 
estimate the line contribution to the filter flux to be ~8% and should not affect the observed 
morphology of the dust continuum. 



3. Results 

3.1. Lobes and Halo 

The HST Ha image of NGC 6644 reveals several morphological structures in the nebula 
(Fig. [1]). Easily discernible are two bipolar lobes (labeled as a — a' and b — b') along the 
approximate SE-NW direction. A ring structure (labeled as "ionized torus" ) near the center 
can also be seen. Beyond the bright ring, there are associated nebulosities (marked as "c"), 
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in particular in the western direction. It is possible that these nebulosities represent the 
projection of a third pair of lobes which are aligned nearly along the line of sight. An 
extended diffuse structure (marked as "halo") can also be seen. In Fig. HI we display the 
image with the central part of the nebula saturated, allowing the fainter outer structures 
to be seen. There are apparently sub-structures in the nebula. In Fig. |5l we mark a few 
possible sub-structures (which appear predominantly in the north) as d,e,f,g,h. Closer to 
the core of NGC 6644 appears several more sub-substructures, again most obvious in the 
northern part of the image. These are labeled i,j, k in Fig. El The axes of these structures 
are defined by the caps of the lobes, location are marked by horizontal bars in Figs. [5]and[6j 
The orientation of the axis of d is less well defined as the cap is not complete, but the others 
are reasonably well defined. What is interesting is that the axis j seems to lie along the same 
orientation as d , as well as the axes of k and g. The fact that these sub-structures appear 
stronger and more evident in the northern half of the image may provide further evidence 
that the north is tilted towards us and therefore experiences less extinction than the south. 

The symmetry axes of lobes a and b intersect approximately at the position of the central 
star. The two pairs of butterfly-shaped, closed-end bipolar lobes (a — a' and b — b') have 
symmetry axes that lie at the position angles (PAs) of -43 ±2° and -27±2°, respectively. The 
lobes a — a' and b — b' have similar projected sizes on the sky. Lobe a — a 1 has a projected 
length of ~ 7. 76" and b — b' has a projected size of ~ 8.81". If we assume a distance of 



2.5 kpc (ICahn et al.lll992l ) the physical size of the total extent of the b — b' lobe is ~ 0.35 
sec# light yr, where 9 is the inclination angle. For an expansion velocity of 100 km s -1 , the 
kinematic age of this lobe is ~ 521 yrs if the lobes lie close to the plane of the sky. From 
this kinematic age, NGC 6644 can be considered as a young PN. 

For the central ring ("ionized torus"), its major axis is found to have an angular orien- 
tation of PA ~ -55±3°. Similar structures have been fou nd in other m ultipolar PNs such as 



the two bright partial rings seen in the waist of He 2-47 (jSahaill2000bl . Fig. lb). The eastern 
segment of the ring in NGC 6644 is brighter than the western side, as evident in both the 
Ha image (Fig. [1]) and in the [S iv] image in Fig. [81 This suggests that the eastern segment 
is closer to us, whereas the western side suffers from extinction from intervening dust in the 
dust torus. 

This bright multipolar nebula is surrounded by a faint, extensive halo with a well- 
defined spherical shape. The extended halo emission at Ha can be seen in an averaged 
radial surface brightness profile (with the two bipolar lobes excluded) shown in Fig. [71 This 
profile is constructed from measurements of 5° intervals between PA=25°-85°and PA=205°- 
265°, after the removal of all field stars. The intensity of the nebula within l."4 is dominated 
by the ionized torus and the angular radius of the faint halo is estimated to be larger than 
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4". The surface brightness of the extended part at radius 4" relative to that in the peak is ~ 
10~ 3 . At a distance of 2.5 kpc, the physical size of the halo is 0.15 light-yr, with a kinematic 
age of ~ 3x 10 3 yr, if the expansion velocity of the halo is 15 km s _1 . 

An intensity profile of the nebula (minus the lobes) averaged over all angles is shown 
in Fig. [7J Also shown is a fit of the profile by a halo density distribution of r~ a , where r 
is the radial distance from the central star, with an inner and outer radii of -Ri n ,out — 0.585 
and 6.5 arc sec respectively. For the modeling, we have assumed that Ha is optically thin 
and its intensity integrated along each line of sight is proportional to the emission measure, 
EM(p) = J n 2 (r)d£, where n(r) is the density of hydrogen ion, p is the angular distance 
from the central star, and di is the path length element along the line of sight. Therefore, 
under the assumption that the halo is a spherical shell, we have 

EM(p) = 2n 2 (R in )R^p 1 - 2a / (cos9) 2a ~ 2 d9, (1) 

Je 1 

where 9i = cos _1 (p/-Ri n ) and for R in < p < R ut and p < i?i n respectively, and 9 2 = 
cos~ l (p/ R ont )- For a stellar wind with a constant mass-loss rate and expansion velocity, a 
has a value of 2 and the density distribution varies as n(r) = n(R m )(r / R\ n )~ 2 . While the 
model provides a reasonably good fit to the halo, it fails in the inner region. 

In a second model, we introduce a face-on torus with a cross-section radius of R t = 0'.'4 
and a homogeneous density distribution of n%. For the torus, we have 

EM(p) = 2n 2 t ^R 2 -(p-R c ) 2 , (2) 

where R c = 0'/585 is the distance from the circle center of the torus to the central star. 
The torus is embedded in a halo with a inner radius of 1 arc sec. This addition makes an 
improved fit to the intensity profile in the central region. The remaining excess in emission 
between 1 and 2 arc sec (as seen in Fig. [7J could be due to emission from the nearly pole-on 
lobes. 

In the model, we do not consider the effect of extinction. Due to the presence of dust 
torus (see below), the central region may have a higher extinction than the halo. If this is 
the case, EM{p) given in eq. (2) could be underestimated and the actual intensity of the 
inner region would be even higher. This will further strengthen the need for a face-on torus. 



-7- 



3.2. The dust torus 

In Fig. M we show the 4-band T-ReCS (PAH-2, Si-5, [S iv], and Q a ) flux-calibrated 
images in units of Jy per square arcsecond. The T-ReCS images of NGC 6644 were each 
deconvolved using the observations of the standard star HD 175775 taken immediately after 
the science observations. The "lucy" task in the "analysis" group of the "stsci" external 
IRAF package was used for the deconvolutions. For each filter a background-subtracted, 
normalized PSF image of size 51 by 51 pixels was created from the standard star observation, 
and this was used as the kernel for the deconvolution. The deconvolution was iterated 25 
times to produce the final deconvolved images. It was found that most of the improvement 
in the images took place within the first 5 to 10 iterations. Under the seeing conditions 
reported in section I2.2[ the Lucy deconvolution is robust to small changes in the FWHM of 
the PSF. The deconvolved images are shown in Fig. [9] together with the Ha image on the 
same scale. 

From Figs. [S] and [HI we can clearly see that there is a spatially resolved structure in 
the central region of NGC 6644. The [S iv] image shows a partial ring structure similar to 
that seen in the central region in the Ha image, suggesting that they both trace the ionized 
gas distribution. The Si5 and PAH2 images are different from the [S iv] image, showing 
a double-peaked structure. The Q a image shows a more diffuse structure than the other 3 
images. The overall maximum angular extent of the emission regions in the four images are 
similar. A more detailed comparison between the HST Ha image with the Gemini [S iv] 
image shows that the ring in the [S iv] image is less complete than its counterpart in the Ha 
image (Fig. [9]). There are two factors that can cause a difference between the two images. 
The [S iv] image includes contributions from both the [S iv] line and dust continuum, and 
therefore should reflect the morphology of the ionized torus as well as the oblique dust torus. 
Since the dust torus (see Fig. [8]) is brighter in the eastern side, the [S iv] image should 
therefore appear as a less complete ring than in the Ha image. If there is also effect of local 
extinction that causes the eastern side to be brighter (see section |3~T|) . then the [S iv] image 
should suffer from less extinction and therefore should show a more complete ring. This 
effect is apparently not as important. 

If the double-peaked structure in the PAH2 and Si5 filter images is interpreted as an 
oblique torus, the size of the torus has dimensions of 28 x 15 pixels with the major axis 
lying along PA ~ -44±3°. This translates to a major axis size of 2'.'52 and a minor axis size 
of 1'.'35. Assuming that this ellipse is a projection of a circle on the sky, an angle of tilt of 
the dust torus is approximately 32° (with 90° being pole on). The fact that this torus is not 
obvious in the optical image is probably because the dust torus is optically thick to the UV 
photons and the volume of the torus is not ionized. 
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In Table 3, we compare the fluxes observed in the four band T-ReCS observations and 
fluxes from the IRS observations. Columns 1 and 2 list the observed filter names and central 
wavelengths of these filters. The corresponding band width for each filter is given in column 
3. The measured integrated fluxes from the T-ReCS observation are given in column 4. The 
simulated T-ReCS in-band fluxes derived from the IRS spectrum and the filter profiles of the 
4 bands are given in column 5. There is general consistency between the Gemini and IRS 
measurements with the discrepancies being due to varying sky conditions and uncertainty 
of photometric flux calibrations. These uncertainties are higher in the Q a band. 



3.3. Infrared features 



The combined IRS spectrum of NGC 6644 shown in (Fig. [3]) shows a strong infrared 
continuum due to dust emission peaking at ~ 30 /im. On top of this continuum are a 
number of emission lines typical of the spectrum of PNs. Among the strongest emission 
lines detected are the fine structure lines of [S IV] at 10.51 /im, [O IV] at 25.88 /im, [Ne III] 
at 15.55 and 36.01 /im, and weaker H recombination lines H I (7-6) at 12.37 /im and H I 
(8-7) at 19.07 /im (Fig. [TO]) . The measured emission line fluxes in the spectrum are given in 
Table 2. The first two columns in Table 2 contain the central wavelengths of the emissions 
and line identifications, respectively. Column 3 gives the observed fluxes measured using the 
Gaussian fitting routine. Using the fine-structure lines [S ill] and [Ne ill], we can derive the 
electron density. From the [S ill] 18.7 to 33.4 /im and [Ne ill] 15.6 /im to 36 /im line 
ratios, we derive logn e = 3.8 and 3.4 cm -3 respectively. The last column gives the total 
flux of the entire nebula by scaling the observed flux by a factor, derived from the average 
theoretical flux ratio between infrared H I lines (H I (7-6), H I (8-7)) and H/3 emission based 
on Case B recombination line theory. Assum ing an electron temp erature of T e = 12,000 K 
and an electron density of n e = 7,200 cm -3 (IShaw fc Kalerlll989l ). the average theoretical 
ratios are then combined to derive the scaling factor of 1.01 for the IRS spectrum. 

Also present is the strong, broad AIB at 11.3 /im. This feature is due to the C-H out-of- 
plane bending of aromatic compoun ds and is commonly seen in carbon-rich PNs and proto- 
PNs (IKwok. Volk. fc Hrivnakill999l ). Several other broad features at 12.6, 14.2, 16.5 /im are 
also likely to be AIB features. A strong, broad 12 /im plateau emission feature, attributed to 
out-of-plane bending modes of a mixture of aliphatic groups attached to the aromatic rings 
(IKwok. Volk. fc Bernathll200ll ). can also be seen. A recent study of AIB features in Galactic 
Bulge PNs using IRS data was presented by lPerea-Calderon et al.l (120091 ). in which they find 
AIB features commonly present in their sample. 
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3.4. Spectral Energy Distribution 



Although PNs were first known for their optical emission-line characteristics, it was 
realized after the IRAS mission that a significant amount of the total energy output is emitted 
in the infrared region due to thermal emission by dust grains. A systematic investigation 
for spectral energy distributions (SEDs) of youn g PNs covering the wav elength range from 
ultraviolet (UV) to far-infrared was performed by lZhang fc Kwokl (119911 ). Depending on the 
stage of evolution, it is found that the photospheric, ionized gas, and dust components are the 
main contributors to the total observed flux. In constructing the SED of NGC 6644, we have 
made use of data in the astronomical data archives of both ground-based and space-based 
observations. 

In the UV region, we have used observations from the International Ultraviolet Explorer 
(IUE) low dispersion spectrograph. All IUE spectra of the object are produced at the 
IUE Data Analysis Center (IUEDAC). Reduced data of this nebula are taken in co-added 
positions of the short wavelength prime (SWP) and long wavelength prime (LWP) spectra 
to enhance the signal-to-noise ratio, respectively. The parts of this spectrum in SWP and 
LWP are then combined at 1975 A to produce a single, resultant spectrum. 

A journal of the IUE data and Spitzer IRS AORkey of the nebula are given in Table 4. 
In the infrared, the Spitzer IRS observations are used together with the photometric observa- 
tions of Infrared Astronomical Satellite (IRAS). The photometric measure ments of NGC 6644 
at the four wavelength bands of 12, 25, 60, and 100 /xm are taken from iTajitsu &: Tamura 
jl998h . 



The B and V magnitudes of the central star of NGC 6644 are observed by lShaw fc Kaler 

( 119891 ) . Near- infrared magnitudes of I, J, H, and K/Ks observed by the Deep Near- Infrared 



Southern Sky Survey (DE NIS) and the Two Micron All S ky Survey (2MASS) are derived 
from DENIS database and iRamos-Larios fc Phillips! (120051 ) . A summary of these archival 
data is given in Table 5. 

After correcting the UV spectrum and optical measurements by an extinction value of 
0.41 (IShaw &: Kalerlll989l ) (which accounts for both circumstellar and interstellar extinction), 
we fitted the emerging flux (F\(total)) by a three-co mponent model includ ing a hot central 
star, a warm gaseous nebula, and a cool dust shell (jZhang & Kwokl Il99ll ). The hot stellar 
component is assumed to be a blackbody of temperature T*. The gaseous nebular continuum 
is the sum of the free-free (f - f), bound- f ree (b - f), and two-photon emissions using the 
emission coefficients given in iKwokl ( 120071 . p. 171) at an electron temperature T e = 12 000 
K and an electron density n e = 7200 cm -3 ( IShaw fc Kalerl I1989I ). The emerging flux is 
therefore given by: 
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F x (total) = F x {s) + F x (g) + F x (d) (3) 

where F x (s), F x (g), and F x (d) are the flux densities of the photospheric continuum of the 
central star, the nebular continuum, and the dust thermal emission continuum, respectively. 
The stellar flux is then given by 



F x (s) = (n%)B x (T,) (4) 

where T* are respectively the angular radii and the effective temperatures of the central 
star, and B X (T*) is the Planck functions for temperatures T*. Given the uncertain extinction 
corrections for the visible photometry of the central star, and the fact that the B and V 
photometry lie in the Rayleigh- Jeans side of the blackbody, the value of T* cannot be deter- 
mined precisely. Our best estimate for is 119 000 ±8500 K, which is con sistent with the 
previ ously reported values of 106 000 K ( Shaw fc Kaler 1989h and 115 000 K ( Zhang &: Kwok 
199ll ). With an assumed distance of 2.5 kpc (jCahn et al.l Il992h . the der i ved lu minosity of 



the central star is ~ 4600 L Q . According to the model of ISchonbernerl (119811 ). the above 



temperature and luminosity of NGC 6644 imply an evolutionary age >6000 years. The value 
is at least two times higher than the observed kinematic age (section l3.ip . This may suggest 
that the visible nebulosity (from which we estimate the kinematic age) originates from very 
late stellar winds, and observations with higher sensitivity are required to detect the fainter 
extended structures. 

The near-IR photometry points suggest that there is excess emission between 2 and 5 
pm which could be due to a cool companion, a hot dust component, or scattering from an 
unseen disk. We have therefore added a blackbody of temperature 2950 K to the model fit. 
Using the same assumed distance, this companion would have a luminosity of 14 L . 

We find that the observed dust emission component is too broad to be fitted by a single 
blackbody. We therefore fit the dust component by blackbodies of two different temperatures, 
a warm (T wd ) and a cold (T cd ) dust components. The total dust thermal emission is therefore 
given by the sum of two terms: 



F x (d) = F x (wd) + F x {cd) 



3M wd Q x B x (T wd ) t 3M cd Q x B x (T ( 



+ 



cd) 



(5) 



Aap s D 2 Aap s D 2 
where M wd and M cd are the masses of warm and cold dust components, a is the grain radius 
which depends on its physical shape, p s is the density of grain, Q x =Q o (X/\o)~ a is the grain 
emissivity function, and D is the distance to the nebula (IKwoki 120071 . p. 312). According to 
our fit, the temperatures of the two dust components are 292 K and 123 K. Assuming Q 



[1 pm)=0.1, p s =l g cm 



=0.1 /im, a =1, and D=2.5 kpc, we obtain 8 x 10 8 M and 
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3.12 x 10~ 5 M for the mass of the warm and cold dust components respectively. Without 
spectral coverage between 3 to 10 /im, it is difficult to constrain the temperature of the warm 
dust component precisely and the above mass estimate for the warm component is therefore 
subject to a large margin of error. 

The SED of the nebula corrected for extinction is shown in Fig. [11] We can see that the 
model gives a reasonable fit to the observed data from UV through the far-infrared. From 
the total observed nebular flux of 5.3 x 10~ 10 erg cm -2 s -1 , an emission measure of 2.5 x 
10 59 cm -3 was derived. At a distance of 2.5 kpc and an electron density of 7200 cm -3 , the 
total ionized gas mass is derived as 0.03 M . From the model fitting, the fraction of total 
fluxes between 0.1 and 100 /zm from the central stellar, nebular gaseous continuum emission, 
and dust components are 38%, 18%, and 44%, respectively. The stellar flux shortward of 0.1 
/im is assumed to have been absorbed by the gas component. 

From the SED, we can see that there is a broad 30 /im feature in NGC 6644. Thi s 
unidentified emis sion feature is commonly seen in carbon stars (yolk- Xiong. fc Kwokll2000r) . 



post-AGB stars ( jSzczerba et al.l Il999t iHrivnak. Volk. fc Kwo 
young PNs. 



2000 



Volk et all 120021 ) and 



4. Simulated 3-D structure of NGC 6644 



The observations of NGC 6644 as reported in section 13.11 suggest that NGC 6644 is 
a new memb er in the club of multipolar nebulae . The clas s of multipolar nebulae , first 
described by iManchado. Stanghellini. fc Guerrero! (119961 ) andlSahai fc Traugerl ( 119981 ) , now 



include many promi nent membe rs including NGC 2440 (ILopez et al. 



19981 1. He 2-47 and 



Ml-3 7 flSahaill200qbh. H e 2-113 (ISahai et al.ll2000d ). NGC 6881 (IKwok fc Sul 120051 ). NGC 
6072 (IKwok et al.ll2010l ). etc. Multipolar nebulae, having a point-symmetric rather than 
axial symmetric structure, need 3-D modelling to visualize their structure. With this aim, 
we have constructed a schematic model of NGC 6644 assuming th ree pairs of bipolar lob es 
and an equatorial dust torus using the software program SHAPE (jSteffen fc Lopez! 120061 ). 



SHAPE is a morpho-kinematic modeling tool intended for the analysis of the 3D geom- 
etry and kinematic structure of gaseous nebulae. SHAPE uses both kinematic and spatial 
observations as a guide for the interactive reproduction of model nebulae. Once a 3D model 
has been constructed, parameters such as position angle and inclination can be manipulated 
in order to study the object from different orientations. In this work we are interested in 
studying the complex 3D morphology of multi-polar PN for which SHAPE is an ideal tool. 
We do not attempt to determine physical parameters such as the temperature and density 
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structure and thus we do not perform any radiative transfer or hydrodynamic simulations. 
Instead, the brightness and location of each component is assigned based on the qualitative 
analysis of the observed images. 

In addition to the observed double bipolar lobes a and b, we hypothesize an additional 
pair of lobes that is perpendicular to the observed ionized torus in Fig. [TJ This pair of 
lobes will be referred to as lobe c. We further assume that all three pairs of lobes are equal 
in length and that the b lobes are oriented in the plane of the sky (inclination 0°). With 
these assumptions, we can then derive an inclination of 31° and -85° for the a and c lobes 
respectively. The estimated symmetry axes of these three lobes are along PA= -42°, -28° 
and 76° for the a, b and c lobes respectively. 

We created a dense torus as the counterpart of the infrared torus seen in the Gemini Si5 
and PAH2 images with a center that coincides with that of the three lobes. If we assume that 
the observed elliptical torus is a tilted circular disk, then the model's major-to-minor axis 
ratio implies an inclination angle of 35° (with 90° being pole on). This value is in reasonable 
agreement with our previous estimate of 32° based on the ellipticity of the torus (see § 3.2). 
In this model, the symmetry axis of lobe a — a' is perpendicular to the infrared dust torus. 
A counterpart of the ring seen in the HST Ha image and the Gemini [S iv] image is also 
included. The model parameters of these structures are listed in Table 6 and a schematic 
mesh model is shown in Fig. | 



This simple model can provide insight into what multipolar PNs look like at different 
orientations. In Fig. [13] we present a three-dimensional representation of NGC 6644 viewed 
from different angles. The bottom left image is the model as seen from Earth, and the 
orientations of the other images are as labeled. The x and y axes indicated are those of the 
image plane and the values represent rotations around those axes in degrees. 

When y = 90°, the nebula reveals three distinct pairs of lobes with an elliptical disk . 



This form is similar to that found in two young PNs, He 2-47 and M 1-37 (jSahail l2000bl ). 



When y = 0° the bipolar lobes appear closer togethe r and overlap to v arying extents. This 



type of structure is se en, for example, in NGC 6072 (IKwok et al.ll2010l ) and the Frosty Leo 



Nebula (jSahail l2000al ). The present model illustrates the possible complex 3-D structures 
that may be present in other objects and indeed it would be difficult to visualize the structure 
of NGC 6644 without the aid of such a model. The main purpose of this simulation is to show 
how complex the structure of a multipolar PN can be and not as a quantitative fit to the 
observations. A proper model of a complex object such as NGC 6644 will require kinematic 
information from all of the morphological features, which can be provided by integral field 
spectroscopy with large optical telescopes. 
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5. Conclusions 

Although NGC 6644 appears to be a typical elliptical nebula, high angular resolution 
and dynamical range optical observations reveal that it in fact has an extremely complicated 
structure. From HST and Gemini observations, we have identified at least two pairs of 
bipolar lobes and an infrared dust torus. The orientation of the dust torus is approximately 
perpendicular to one pair of the lobes. In addition, there is a ring of ionized gas lying almost 
on the plane of the sky. We suggest that there could be a 3rd pair of bipolar lobes whose 
axis is perpendicular to this ring. Several other faint sublobes at different radial distances 
and angular directions can also be seen. 

It is clear that NGC 6644 is a multipolar nebula. In order to illustrate its possible 3-D 
structure, we have constructed a model to show what it may look like when viewed from 
different perspectives. Its simulated rotated images resemble some other multi-polar nebulae, 
suggesting that nebulae with different appearances may have similar intrinsic structures. 

The emergence of the class of multipolar nebulae has greatly altered our perception of 
the morphological structures of PNs. The commonly assumed simple structures of PNs are 
probably the result of inadequate sensitivity and spatial resolution imaging. Deep imaging 
of PNs is needed to reveal the true intrinsic structures of PNs. 
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RA [2000] 

Fig. 1. — HSTHa image of NGC 6644 showing various morphological structures. The inner 
(central) and outer regions are displayed at different levels to better show the structures. 
The two bipolar lobes are marked as a — a' and b — b' and the central ring structure is marked 
as "ionized torus". The nebulosity associated with the ring (marked as "c") could be the 
projection of another pair of lobes aligned almost along the line of sight. An extensive halo 
(marked) can also be seen. 
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Fig. 2. — Apertures of the Spitzer IRS observations overlaid on the HST image of NGC 6644. 
The red, blue, and green boxes show the SH (10-20 jum), LH (19-37 /jm) and SL (5.2-14.5 
/jiii) observations respectively. 
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Fig. 3. — Comparison between the Spitzer IRS spectra and the Gemini T-ReCS photometry. 
The observations of the IRS SH and SL have been adjusted upward by factors of 1.29 and 
1.625 respectively (see text). The T-ReCS and IRAS photometric measurements are shown 
as circles and squares respectively. The Gemini T-ReCS filter transmission functions are 
plotted as dotted lines. 
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Fig. 4. — The HST Ha image of NGC 6644 with an intensity display setting to show the 
outer structures of the nebula. The intensity display is on a logarithmic scale and the grey 
scale bar is given at the bottom in units of counts per pixel. 




Fig. 5. — A number of faint sub-lobes can be seen in the HST Ha image of NGC 6644. The 
caps of the lobes have been highlighted by bars for easy identification. The lines perpendic- 
ular to the bars labeled as d, e, f, g, h are inferred axes of the lobes. 
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Fig. 6. — Three inner lobes are marked on the HST Ha image of NGC 6644. The labels are 
as in Fig. 
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Fig. 7. — Profile of surface brightness distribution of Ha emission in NGC 6644 averaged 
over all angles except the regions occupied by the lobes a — a' and b — b'. The halo can be 
fitted by a density law of r -2 (model in solid line) but the intensities in the inner regions 
require contribution from a separate component (dashed line). See text for details. 
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Fig. 8. — Flux-calibrated Gemini T-ReCS images of NGC 6644. The four panels give 
respectively the PAH-2 11.3 jum image (upper left), the Si5 11.66 /zm image (upper right), 
[S IV] 10.52 /im narrow-band image (lower left), and Q a 18.3 /im medium-band image (lower 
right). North is up and East is to the left. The intensity display is on a linear scale and the 
grey scale bar is given at the bottom in units of counts per pixel. 
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Fig. 9. — Comparison between the HST H a image (lower left) and the flux-calibrated, de- 
convolved, false-color Gemini T-ReCS images displayed on the same pixel scale. The Gemini 
[S IV] 10.52 /im narrow-band, PAH-2 11.3 /im narrow-band, and the Q a 18.3 fim medium- 
band images are given in the upper left, upper right, and lower right panels respectively. 
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Fig. 10. — Spitzer IRS spectrum of NGC 6644 in wavelength range from 10 /im to 37 /im. 
The emission lines and AIB features are marked. The Gemini T-ReCS filter transmission 
functions are plotted as dotted lines. 
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Fig. 11. — SED of NGC 6644 in the wavelength range from 0.1 /im to 170 /im. The filled 
squares are B and V photometry of the central star, the filled- circles are DENIS near-IR 
photometry, the open triangles are 2MASS photometry, and open circles are IRAS 12, 25, 
60, and 100 /im photometry. Note that the flux measured from IRAS 100 /im is an upper 
limit. The four BB-like curves (from left to right) represent the central star, cool companion, 
and the two dust components. The model curve with step jumps represents the nebular 
component. The total flux from all components are plotted as a solid line on top. 



Fig. 12. — Comparison of optical-infrared observations and the three-dimensional model of 
NGC 6644. (a) The gray frame shows the HST Ha image, as in Figure 1. Linear contour 
plot of the T-ReCS PAH-2 11.3 fim image was overlaid on the HST image (red lines), (b) 
SHAPE three-dimensional mesh model. The bipolar lobes a, lobe b, and lobe c are displayed 
in red, blue, and green, respectively. The equatorial ring observed in the HST Ha image is 
shown in light blue and the infrared torus observed in Gemini PAH2 image is in white. 
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Fig. 13. — The model of NGC 6644 as viewed from different angles. The rotation angles 
around the x and y axis (of the image plane) are given in the top corner of each image. The 
bottom left image represents the model as viewed from Earth (i.e. no rotation). The color 
scheme is the same as in Fig. [12] 
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Table 1. Summary of T-ReCS Observations of NGC 6644 



Observation Date 


Filter Name 


A c 


A A 


Exposures 






(cm) 


(Mm) 


(<0 


2007 Oct 10 


[SIVJ-10.52 fim 


10.52 


0.17 


455.1 


2007 Oct 10 


PAH2-11.3 fim 


11.30 


0.61 


461.5 


2007 Oct 10 


Si5-11.66 urn 


11.66 


1.13 


238.9 


2007 Oct 10 


Qa-18.30 nm 


18.30 


1.38 


238.9 



Table 2. Measured emission line fluxes in the IRS spectrum 



Wavelength 


Identification 


Observed flux 


Normalized flux 


W 




(10- 10 crgs cm- 2 s- 1 ) 


(10- 10 crgs cm- 2 s- 1 ) 


10.51 


[S IV] 


0.1098 


0.1109 


11.20 


AIB 


0.1513 


0.1530 


11.76 


[CI IV] 


0.0012 


0.0013 


12.37 


HI (7-6) 


0.0024 


0.0025 


12.81 


[Nc II] 


0.0055 


0.0056 


13.09 


[At V] 


0.0019 


0.0020 


14.32 


[Ne V] 


0.0225 


0.0227 


15.55 


[Ne III] 


0.2366 


0.2389 


18.71 


[S III] 


0.0157 


0.0159 


19.07 


HI (8-7) 


0.0016 


0.0017 


24.32 


[Ne V] 


0.0087 


0.0088 


25.88 


[O IV] 


0.1440 


0.1454 


33.47 


[S III] 


0.0038 


0.0039 


36.01 


[Ne III] 


0.0198 


0.0200 
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Table 3. Comparison of IRS and Gemini T-ReCS Observations for NGC 6644 



Filter Name 


Central Wavelength A c 


Band Width A A 


Flux Density 


IRS Flux Density" 




(/mi) 


(/im) 


(Jy) 


(Jy) 


SIV-10.52 /urn 


10.52 


0.17 


4.18 


4.29 


PAH2-11.3 /tin 


11.30 


0.61 


2.31 


2.41 


Si5-11.66 /tm 


11.66 


1.13 


2.39 


2.32 


Qa-18.30 /im 


18.30 


1.38 


5.78 


5.01 



a The IRS flux density is obtained from integrating the convolution of filter transmission curve and Spitzcr 
IRS spectrum. 



Table 4. Available IUE data and IRS observation 



PN 


Name 


Instrument 


Exposures 


IUE Spectra 


G 008.3-07.3 


NGC 6644 


SWP 31711 


1800 






SWP 01734 


1800 






LWR 01630 


2400 


IRS Observation 


G 008.3-07.3 


NGC 6644 


AORkcy 11334400 


347 
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Table 5. Other photometric measurements of NGC 6644 



Filters 


Flux/Flux density 


Central star 


B (mag) a 


16.6 


V (mag) a 


15.63 


-log F (H/3) (ergs.cm^ 2 ^- 1 )" 


11.01 


Nebula 


DENIS I 6 (mag) 


12.678 


DENIS J b (mag) 


11.743 


DENIS K b (mag) 


10.695 


2MASS J c (mag) 


11.562 


2MASS H c (mag) 


11.599 


2MASS Ks c (mag) 


10.627 


IRAS F 12 iim d (Jy) 


2.19 


IRAS F 25 fj,m d (Jy) 


7.57 


IRAS F 60 fj,m d (Jy) 


4.22 


IRAS F 100 nm d (Jy) 


18.39 



"From Shaw & Kaler (1989); 

b From DENIS database; 

c From Ramos-Larios & Phillips (2005); 

d From Tajitsu & Tamura (1998), the color- and 
diameter-corrected IRAS fluxes were given. Note that 
the flux measured from 100 (im is an upper limit. 



Table 6. Comparison of Observed and Model Parameters of the Lobes and Torus 



Value 



Parameters 


Observed 


Model 


Comment 


Position angle of lobe a — a 1 


-43°±2° 


-42° 




Position angle of lobe 6 — 6' 


-27°±2° 


-28° 




Position angle of lobe c — d 




76° 




Inclination angle of lobe a-a' 




31° 


Assuming the orientation angle of sky plane is 0° 


Inclination angle of lobe 6 — 6' 




0° 




Inclination angle of lobe c — c' 




-85° 


11 


Position angle of dust torus 


-44° ±3° 


-41° 


Derived from peak position of torus 


Inclination angle of dust torus 


32° 


35° 


Derived from the major-minor axis ratio 


Position angle of ionized torus 


-55°±3° 


-60° 




Inclination angle of ionized torus 


-106° 


-100° 





